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Abstract. We compare and contrast two different metric based formulations of non- 
linear cosmological perturbation theory: the MW2009 approach in [K. A. Malik and 
D. Wands, Phys. Rept. 475 (2009), 1.] following Bardeen and the recent approach 
of the paper KN2010 [K. Nakamura, Advances in Astronomy 2010 (2010), 576273]. 
We present each formulation separately. In the MW2009 approach, one considers 
the gauge transformations of pcrturbative quantities, choosing a gauge by requiring 
that certain quantities vanish, rendering all other variables gauge invariant. In the 
KN2010 formalism, one decomposes the metric tensor into a gauge variant and gauge 
invariant part from the outset. We compare the two approaches in both the longitudinal 
and uniform curvature gauges. In the longitudinal gauge, we find that Nakamura's 
gauge invariant variables correspond exactly to those in the longitudinal gauge (i.e., 
for scalar perturbations, to the Bardeen potentials), and in the uniform curvature 
gauge we obtain the usual relationship between gauge invariant variables in the flat 
and longitudinal gauge. Thus, we show that these two approaches are equivalent. 
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1. Introduction 

Many problems in Physics and Applied Mathematics can be described by non-linear 
systems of evolution equations. These are notoriously difficult to solve exactly because 
of the non-linearity. An example of such a theory is General Relativity. Einstein's 
equations are highly non-linear and can only be solved exactly in a small number of 
useful cases. To go beyond these solutions perturbative methods are used. Given a 
solution to the equations in the form of a metric gj^J we assume that we can approximate 
a neighbouring, more general, solution g^j^^, using a power series. Thus we express the 
more general solution in the form 

1 

2- 



9,u = gl^J+gl!J + ,:g^^J + ... . (i-i) 



9,u = g'^^J+egl!J + -e'g^^J + ... , (1.4) 



The metric gj^J is called the background and the remaining terms are perturbations of 
the background. The ffist order part is 

9,u-g^^J^gl'J, (1.2) 

where the remaining terms are assumed to be negligible compared to gju! and they 
are neglected at first order. In a similar way the higher-order perturbations can be 
identified, so at second order we have 

gfiu ~ g^u ~ g^iv — g^v ■> 

(1.3) 

and so on. This can be described simply if we introduce an infinitesimal parameter 
e ^ 1 for the perturbation and assume that the series can be written as 

1 

— ( 

2 

where the quantities with bars have absolute magnitudes less than one. In this format 
the orders correspond to the powers of e. In practise it is often a nuisance to introduce the 
parameter e so we will use the form (11. ip if there is no confusion. Issues of convergence 
can, in general, be removed by working in a small enough neighbourhood but they 
should not be ignored. 

Having set up the approximation (II. ip we have to substitute it into the Einstein 
equations 

Gf,u + ^gf,u = Tf,^ , (1-5) 

to obtain solutions of the required order of approximation. This is more complicated 
than might be expected. Perturbations of the metric imply perturbations of the energy 
momentum tensor and vice versa but more significantly, calculation of the connection 
coefficients and the Ricci tensor involve raising and lowering indices and so introduces 
more terms and potentially couple terms of different orders. At zeroth and ffist order 
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this is not a problem but at higher orders it makes calculations much more complicated. 

_(2) 

Even at second order there are "proper" second order terms, for example gl^J , and terms 
quadratic in the first order quantities, for example gj^J . 

In addition to the complexity of the system of equations, the split of the metric and 
matter variables into a background and perturbations introduces spurious coordinate 
artefacts or gauge modes as described in detail in numerous reviews [H [21 [3l HJ [5] . We 
only give a brief explanation, here. 

As we are interested in cosmology, we will assume that the real universe is described 
by a perturbed Friedmann-Robertson- Walker (FRW) metric and that the background, 
or unperturbed spacetime, is described by an exact FRW metric. So we have in effect 
two spacetimes - one "physical" and one "fictitious" . We label points in the background 
by coordinates {x^} and a one-to-one- map between points in the background and points 
in the physical spacetime maps these coordinates from the background to the physical 
spacetime. We refer to this one-to-one-map as a gauge choice. A change in the map is 
called a gauge transformation and this may be carried out in a number of ways, see for 
instance the recent reviews [HI IH E]. A perturbation of some quantity is the difference 
between the value at a point in the physical spacetime and the value at the point in 
the background with the same coordinates. Clearly such a perturbation depends on the 
above gauge choice. 

It is important to note that a gauge transformation is different to a coordinate 
transformation which changes the labels on points in the physical and background 
spacetimes together, and so it does not change the gauge. A simple example of a gauge 
transformation is an implementation of a coordinate change in the physical spacetime 
but not in the background. This changes the correspondence between the points in the 
two spacetimes, so it is a gauge transformation. It is easy to see that a scalar, e.g. the 
energy density p, which is (at least) time dependent, will not be invariant under such a 
transformation. Furthermore if a gauge is chosen to simplify the metric on the physical 
spacetime and some residual gauge freedom remains then spurious gauge mode solutions 
may appear. For these reasons gauge invariant formulations of perturbations and other 
special gauge fixing forms have been sought. These fall into two broad classes, (1) those 
following the general pattern of Bardeen [3 [H El |3l IH |5] and (2) those following the 
covariant form developed by Ellis and Bruni [HI |9l [TOl [TTj . 

At first order in the perturbations the approaches in class (1) above differ largely due 
to different splits of the spatial part of the metric, and notation and sign conventions. At 
second order things are more complicated and there are at least two different approaches. 
Relating the approach used by Malik and Wands |1] and the Nakamura approach 
[5] (hereafter referred to as MW2009 and KN2010, respectively) is the aim of this 
paper. We aim at highlighting similarities and differences of the two subclasses of 
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approaches following Bardeen and try to keep the mathematical background discussions 
to a minimum, referring the interested reader to the reviews listed above and the original 
literature cited therein. However, in order to introduce some of the quantities and 
concepts used later on, we here briefly review the formulation of perturbation theory in 
a more rigorous sense. 

In relativistic perturbation theory, we consider two distinct spacetimes: the 
'background' spacetime, denoted A4o and the 'physical' spacetime denoted ^A, following 
Stewart [26]. The physical manifold is nature itself, and we want to describe the 
properties of this spacetime through the perturbations. On the other hand, the 
background spacetime is a mere reference spacetime on which to calculate perturbations. 
We then introduce a point identification map X : M.q M. relating points on the 
unperturbed manifold to points on the perturbed manifold. By virtue of the pullback of 
this point identification map (denoted A**), we may treat a tensor field on the physical 
spacetime T as a tensor field on the background spacetime A'*T, which we will often 
denote, equivalently, as xT^- The choice of map X, the gauge choice, is not unique; 
if we choose another different map, 3^, the pulled back variables are then different 
representations of the tensor T. These two gauge choices induce a diffeomorphism 
$ : Mq Mq, defined as $ := X'^ o y. 

We conclude this introduction by presenting some essential equations from the 
reviews [HI H], in the notation which is used in the standard literature (e.g. Refs. [H El 
[3]). As in the case of the metric tensor Eq. fll.ip above, we assume that any tensorial 
quantity T can be expanded into a power series by 



where the subscripts denote the order of the perturbation. The change in a perturbed 
quantity (to a certain order) induced by a gauge transformation is given by the 
exponential map, once the generating vector of the gauge transformation, has been 
specified. The exponential map is 



where £^ denotes the Lie derivative with respect to Expanding the exponential map 
and using Eq. (II. 6p we obtain in the background, at first order and at second order, 
respectively. 



T =To + 5T, 

5T = Ti + + ^Tg + . . . , 



(1.6) 



<I>*T = T = e^fT 



(1.7) 



— 



Ti — Ti + f^^To , 



(1.8) 
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where ^'^ is the vector field generating the transformation and is expanded order by 
order as ^'^ = e^f + |e^^2 + ^(^'^)- The exponential map can also be applied to the 
coordinates x'* to obtain the following relationship between coordinates at two points, 
p and q 

x^'{q) = e^^^lp . (1.9) 

Expanding to second order gives 

x'{q) = x^(p) + eer(p) + le'[^UP)^i(P)+^2ip)) ■ (1-10) 

Note that this coordinate relationship is not required to perform calculations in the 
active approach, but will be useful for later discussion. 

On the other hand, when evaluating the gauge transformation rule, a different 
point of view is adopted in the formulation in KN2010, where the change in a perturbed 
quantity induced by a gauge transformation $ is represented in the general form of the 
Taylor expansion of $* 

T(g) = ($*T)(p) = Tip) + e £^,T\^ + ^e' {£^, + 4) T|^ + 0(e3).(l.ll) 

As shown by Bruni and coworkers fTZl [T^ [TU [T51 [IS], the Taylor expansion of the 
pull-back of tensor field is always given in the form of Eq. (11. 6p . even if $* is not an 
exponential map. In this sense, the Taylor expansion (11.111) represents the perturbative 
expansion of a wider class of diffeomorphisms than exponential maps. Through this 
general formula (II. lip and the perturbative expansion (11.61) of the variable T, we can 
reach the same order-by-order gauge transformation rules as Eq. (II. 8p . 

Although the gauge transformation rules at each order appear to have the same 
form in the two different formulations [HE], we should point out that the approaches to 
obtain the transformations are conceptually different. As mentioned above, the Taylor 
expansion (II. lip is valid for a wider class of diffeomorphisms beyond the exponential 
map. Although we may regard Eq. (II. lip as that of an exponential map through a special 
choice of (s-g- Ref. [I2]), there is no guarantee that this is true for any choice of ^2- 
However it seems that, when working to a particular order in perturbation theory, one 
can always use an exponential map that generates the specified gauge transformation 
to that order. Therefore, while conceptually the gauge transformation rules used by 
the different formulations are not the same, this is only a philosophical issue and, in 
practice, the two are equivalent. 

In this paper, we discuss the equivalence of these two formulations in MW2009 
and KN2010, through clarifying the correspondence of the variables. Having introduced 
the basics, we now go on to tackle the main topic of this article, namely, obtaining 
the relationship between the second-order metric perturbation in MW2009 and that in 
KN2010. It should be noted that the two formulations use very different notation. We 
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should note that they use very different notation in these two formulations. Rather than 
try and enforce a common notation between the two, we present each in its conventional 
notation, relating them to one another at the end. Further, we discuss two different 
gauge fixing. One is the Poisson gauge (longitudinal gauge) choice and the other is the 
fiat gauge choice. In MW2009, they showed that these two gauge fixings are complete 
gauge-fixing, while there is no explicit gauge- fixing in the formulation in KN2010. We 
clarify the correspondence between the variables in MW2009 and KN2010 through these 
two gauge-fixing. 

We also have to emphasise that although this paper is not a complete survey of 
the many different formulations of second-order cosmological perturbation theory, the 
correspondence which is clarified here is a useful check of the equivalence of different 
formulations and will hopefully lead to a consensus in the community. 

This paper is organised as follows. In the next section, we describe first order 
perturbations, first in the MW2009 approach and then in the KN2010 approach. In §3] 
we review second order perturbation theory, again starting with the MW2009 approach 
and then the KN2010 approach. In ^ we compare the two approaches, in both the 
longitudinal (or Poisson) gauge, and the uniform curvature gauge. Finally, we summarise 
our results in §3 

If not otherwise stated we use conformal time, t], related to coordinate time t by 
dt = adrj, where a(?7) is the scale factor, see Eq. (12. ip . throughout. Derivatives with 
respect to conformal time are denoted by a prime, and the Hubble parameter is defined, 
in terms of conformal time, a.s H = a' /a. Greek indices, u, A, run from ... 3, and 
lower case Latin indices, i,j, k, run from 1 ... 3. We also make use of abstract indices 
a, b, c, in parts. 

2. First order cosmological perturbations 

The background spacetime Mq considered in cosmological perturbation theory is 
a homogeneous, isotropic Friedmann-Robertson- Walker universe foliated by three 
dimensional hypersurfaces S(?7), parametrised by conformal time r]. In this paper we 
restrict ourselves to considering fiat spatial hypersurfaces, and the line element for this 
spacetime is 



where a = a{r]) is the scale factor and 6ij is the metric on the fiat space. The full 
spacetime metric is then expanded, as in Eq. (II. 6p . as 



ds'^ = a^{ri) — drj^ -|- dijdx^dx^ , 




9^u = 9^ 
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in the notation of MW2009. Alternatively, one can represent the background spacetime 
metric as 

gab = 0^(11) - {dr])a{dr])b + Sij{dx')a{dx^)b , (2.2) 
in the abstract index notation with the full spacetime metric then being expanded as 

gab = g^ab + ^Kb + \eXb + 0(e3). (2.3) 

Equations (I2.2p and fl2.3p are the abstract index notation in KN2010. 

The metric components at each order in perturbation theory can then be expanded 
into scalar, vector and tensor components, according to their transformation behaviour 
on spatial hyper surf aces. 

2.1. MW2009 formulation 

At linear, or first, order in perturbation theory the general scalar, vector and tensor 
perturbations to the flat [K = 0) FRW background spacetime can be expressed in the 
line element 

ds"^ = a^{ri) [-{1 + 2(pi)dr]'^ + 2Budx'dT] + (5^ + 2Cuj)dx'dx^] . (2.4) 

The perturbations of the spatial components of the metric can then be further 
decomposed as |1] 

Bii = Bi^i — Sii , (2.5) 

Ciij = -tpiSij + Ei^ij + Fi(ij) + ^hiij , (2.6) 

where (f)i,Bi,tlJi and Ei are scalar metric perturbations, Su and Fu are divergence- 
free vector perturbations, and hij is a transverse, traceless tensor perturbation. In the 
notation of KN2010, this is then 

h^^ = -2a^<Pi, (2.7) 

hir, = a^Bu = a^DiBi - a^Su, (2.8) 

hij = 2a^Cuj = 2a^ (^-^JiS^j + DiD^Ei + + ^hu,^ , (2.9) 

where Di is formally the covariant derivative associated with the spatial metric 6ij 
and, in practice for this work, it reduces to a partial derivative denoted by a comma. 
Consistently perturbing the spacetime will naturally invoke perturbations to its matter 
content as well. In the following, however, we shall only use the energy density. 

Before studying the transformation behaviour of perturbations at first order, we 
split the generating vector C,i into a scalar temporal part ai and a spatial scalar and 
divergence- free vector part, respectively Pi and 71*, as 

er = (ai,/3i; + 7iO- (2.10) 
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We can then consider transformations of different types of perturbation independently, 
since they decouple at linear order. For example, Eq. f ll.8p implies that the energy 
density perturbation transforms, at first order, as 

6pi = 6pi + p'qOi, (2.11) 

where we have used the fact that the Lie derivative, when acting on a scalar function, 
is just £^ = ^^'•{d / dx^''). The transformation behaviour of the metric tensor, noting that 
the Lie derivative for a type (0,2) tensor is given by 

£i9i.. = 9^.u,xt + g^.xi'^u + ^7A.e,\ , (2.12) 

is 

= Sg^J + CM + . + 9^ • (2-13) 

We can obtain the transformation behaviour of each particular metric function by 
extracting it, in turn, from the above general expression using the method outlined 
in e.g. MW2009. As mentioned above, we do not focus on details here but instead 
quote the results. The scalar metric perturbations transform as 

01 = 01 + Hai + a[ , (2.14) 

^ = ^/>i-i7ai, (2.15) 

= 5i + (2.16) 

^ = ^i + /3i, (2.17) 
the vector metric perturbations as 

^* = 5i*-7/, (2.18) 
Fi^ = Fi* + 7i% (2.19) 

and the tensor perturbation, huj, is gauge invariant. Finally, the scalar shear, which is 
defined as ai = E[ — Bi, transforms as 

ai = ai + ai , (2.20) 

which will be useful later when we come to define gauges in ^ 

2.2. KN20 10 formulation 

An alternative approach within perturbation theory was presented in KN2010 [S], where 
the procedure proposed by KN in 2003 [18] is used to construct gauge invariant variables. 
Evaluating Eq. fll.llj) to first-order perturbation, the gauge transformation rule of the 
first order metric perturbation hab is given by 



hab =■ £^igab- (2.21) 
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We decompose the linear metric perturbation, hat, as 

hab ='■ T^ab + £xgab, (2.22) 

where Hab and £xgab are the gauge invariant and variant parts of the first order 
metric perturbations [H] , respectively and X"' is defined below. That is, under a gauge 
transformation, these are transformed as 

yHab — xHab = 0, yX"" — xX"' = C,^iy (2.23) 

We note that the decomposition in Eq. (12.221) is an assumption, however for 
perturbations to a FRW spacetime it can be shown to be correct 0. 

To proceed, we consider the scalar- vector-tensor decomposition of the components 
hirij hij of hab as 

hi^ = DihiyL) + h(y)i (2.24) 
hij = \^{L)5ij + ^DiDj - h(TL) + 2Diih(Tv)j) + /i(TT)ii|(2.25) 

where A := -D*-Dj = 6^^DiDj. Further h(Tv)j, and h(^TT)ij satisfy the properties 

D'h^v)i = 0, D'h^TV)i = 0, h^T)' , := S^'^t),^ = 0, (2.26) 

h{TT)ij = h(TT)ji, D^h(^TT)ij = 0. 

The generator of the gauge transformation, is also decomposed as 

= Udv)a + + ^iT)^) {dx')a. D'^^T)^ = 0. (2.27) 

Using Eq. (12.211) we can then obtain gauge transformation rules for the components of 

hab- 

yhrir, — xhrjij = 2 (9,, — H) ^rf, (2.28) 

MvL) - xhvL) =^, + (9, - 2H) ^(L), (2.29) 

2^(v)i — xh(v)i = (dr, — 2H) ^(Qn)j, (2.30) 

- a^Aht^L) = - 2H^, + '^A^^, (2.31) 

a^:ih(TL) — Ah{TL) = 2^(L), (2.32) 

o?:^(TV)i — a^Ah{TV)i = ^{T)i, (2.33) 

a^yk(TT)ij - a^:^{TT)ij = 0. (2.34) 

We then inspect these gauge transformation rules and define gauge invariant 
variables. Firstly, Eq. fl2.34p shows that the transverse-traceless part /i(Tr)ij is itself 
gauge invariant, as expected. We denote this as 

(1) 

Xij:= h(TT)ij, (2.35) 
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Secondly, the gauge-transformation rules fl2.30p and ( I2.33P give the transverse vector- 

(1) 

mode Ui defined as 

% := h(^v)i - a^dr^h(Tv)i- (2.36) 

In addition to the vector and tensor modes there are two scalar modes in the first order 
metric perturbation, hab- To see this, we first consider the gauge transformation rules 
fl2.29p and fl2.32p . From these transformation rules, the variable defined by 

:= — -O^dr^hiTL) (2.37) 

transforms as 

yXr^ — x^ri = ^ri- (2.38) 

Using this definition of X^, and inspecting the gauge transformation rule (12 .28 p . we can 

(1) 

show that the variable $ defined as 

- 2a^ $ := hr,^ -2(3^- H) (2.39) 

is gauge invariant. Furthermore, from gauge transformation rules (12.311) . (12.321) . and 

(1) 

(I2.38p . the variable \1/ defined by 

- 2a^ (^^L) - ^A/i(rL)^ + (2.40) 

(1) (1) (1) (1) 

is gauge invariant. The set of variables {$, ^, i/j , Xij} is the complete set of gauge 
invariant variables. 

We can now write the original metric perturbation, hab-, in terms of these gauge 
invariant variables as 

hr,ri = - 2a^ $ +2 (Or, - H) X^, (2.41) 

hr^i = +a^dr^h(TV)i + Dih^vL), (2.42) 

h,j = - 2a^ Sij + Xij WD^D^h^TL) - 2HX^Sij + 2a^D^ih^Tvi3-^^) 

(1) (1) (1) (1) 

From the gauge invariance of the variables Xij, z^j , $, and we may read off the 
gauge invariant part of hab as 



y-ab = a^ {-2 $ {dr])a{dr])b + 2 u\ {dr])^a{dx' 



+ ( -2 5ij+ X^j ) {dx')a{dx% \ , (2.44) 
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The remaining gauge dependent parts in Eqs. fl2.4ip - fl2.43p should then be given in the 
form £xgab for some vector field Xa- In fact, such a vector field is 

Xa := Xr,{dr])a + (^h^TV)i + ^A/i(TL)^ {dx')a, (2.45) 

where X^ is defined in Eq. f l2.37p . One can also check that this vector field Xa satisfies 
Eq. t[223]). 

Thus, it has been shown that the decomposition fl2.22p for the first order metric 
perturbation is correct for cosmological perturbations. We should note that, in order to 
accomplish Eq. f l2.22p . we have assumed the existence of the Green's function A~^, and 
that the perturbative modes which belong to the kernel of the operator A have been 
neglected. However, we ignore both of these issues in this paper, pointing the diligent 
reader to, for example, Refs. [T^ I2U] for more information. 

Finally, the correspondence between the variables for linear perturbations for 
the KN2010 approach presented in this section and those for the MW2009 approach 
presented in § 12.11 are 

hrir, ^ — 2a^0i, h(^vL) ^ a^-Bi, h(v)i —o^Su, 
2 

^ -2^1 + 2 AEi, hiTL) ^ 2Ei, h(Tv)i ^ ^u, h(^TT)ij ^ fiM^) 

Note that at first order there is no difference between the two formalisms, as can be 
seen from Eq. fl2.46p . 

3. Second order cosmological perturbations 

Having reviewed linear cosmological perturbation theory in both the MW2009 approach 
and that in KN2010 above, in this section, we review the formalisms at second order. 

3.1. MW2009 formulation 

In order to consider cosmological perturbation theory to second order we do not truncate 
the perturbative expansion after the first term in Eq. fll.ip . The metric tensor then has 



components form 

55£) = -2a^(r/)</>2, (3.1) 

5g^S = o'{v){B2,-S2d. (3.2) 

= ci\v) ( - 2^25^i + 2^2,^, + 2F2(i,,) + h2i^ , (3.3) 

where the quantities here are analogous to their first order counterparts. At second 

order we split the generating vector ^2 > as at first order, as 

e2^ = (a2,/32; + 72T (3.4) 



Comparing different formulations of non-linear cosmological perturbation theory 



12 



Then, using Eq. (11. Sp . we find that the second order energy density perturbation 
transforms as 

Sp2 = 5p2 + Poa2 + ai(Poai+Poa'i + 25p;) + (25pi + p[,ai),fc(/3i,^+7i^) , (3.5) 

where we note for the first time here that, while at linear order different types of 
perturbation (scalar, vector and tensor) decouple, this is no longer true at higher order. 
This is a crucial qualitative difference between first and second order perturbation theory 
and can lead to the generation of, for example, second order gravitational waves [21] or 
vector modes and vorticity [22]. At second order the metric tensor transforms as 



+ ^?1a (ei + i\ .) + gf^ (ei + ,) • (3.6) 

From this we can extract, as at first order, the transformation behaviour of the metric 
perturbations. Again, we refer to MW2009 [1] for the details and only quote the results 
in this section. One finds that the scalar metric functions transform as 

02 = 02 + Ha2 + a2' + ai [a/' + Si/a/ + {H' + 2H^) ai + 4f/0i + 20'(|.7) 



+ 2ai' (a/ + 20i) + (a/ + Ha^ 



Bo 



2^ k 



52-«2 + /32 + V-Vfc, 



(3.8) 

^VX\ , (3.9) 

(3.10) 



where Xj, and Xbi are defined as 



Xp 



(3.ii; 



= 2 
+ 4 



+ 4/7ai (^i.j + 6j,i) - 2ai,,aij + 2^ifc4/' + ai (C^.^^. + ^^,) 
+ J + iij,kii ,i + + i'ijOii,i + + ^ij,ifc) ^i'' • (3.12) 

We furthermore find that the second order metric vector perturbations transform as 

S2i = S2i - 72/ - Xbj + V^Xe'^tj , (3.13) 



2^kl 



i^2. + 72. + V-%,'=-V-^V-^X 



,kli 5 



(3.14) 
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and the tensor perturbation as 

h2^, = h2^, + X,, + 1 (V-^X'^',,, - X\) 5,, + ^V-^V^^X'^',,,,^. 

+ - (X,,,^ ^ + X^.,,'^ ,) . (3.15) 

Note that the second order tensor perturbation, h2ij, changes under a gauge- 
transformation [B], unhke the first order tensor perturbation, hnj. However, we show 
that the expression (13.151) is itself gauge- invariant as discussed at the end of the next 
section. Hence, using Eq. (I3.15p . we can construct gauge-invariant tensor perturbations 
at second order 



3.2. KN20 10 formulation 

We now consider the construction of the gauge-invariant variables for second-order 
metric perturbation lab in KN2010 [5]. As we have confirmed the decomposition of the 
first-order metric perturbation, hah, into the gauge-invariant l-Lab and gauge- variant parts 
Xa, we can now also find gauge invariant variables for higher-order perturbations |18j . 
Using Eq. (11.81) . the gauge transformation rule for lab is 

ylab — xlab = + '2£^ixhab + {£(2 + <^|i) 9ab- (3.16) 

Through the first order gauge- variant variable Xa, we first define the tensor field Lab as 

xLab '■= xlab — '2£ xXxhab + £l,x9ab (3-17) 

so that it transforms as yLat — xLab = £agab, where a" := ^2 + [^i^A'^]'^- Since a" 
is an arbitary vector field on A4o, this gauge-transformation rule is the same as that 
for the first order metric perturbation, hab (12.211) . Using a similar procedure to that 
with which we decomposed hab, the variable Lab is then Lab ='■ i^ab + £Ygab- Thus, the 
decomposition of the second order metric perturbation, lab, is 

lab ='■ ^ab + '2£xhab + {£y — £x) 9ab, (3.18) 

where Cab and Y°- are the gauge invariant and variant parts of the second order metric 
perturbations, i.e., 

yCab - xCab = 0, yY- - = ^^.^ + [^(1), ^X^. (3.19) 

We may then choose the components of the gauge invariant variables Cab in 
Eq. (I338|) as 

f (2) (2) 

Cab = -2 $ {dr])aidr])b + 2 z/^ {dri)(^a{dx\) 

( (2) (2) \ .1 
+ f -2 * XiA {dx')a{dx')b \ , (3.20) 
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(2) (2) 

where v i and Xij satisfy the equations 

• (2) • (2) 

^^4 = 0, x% = 0, D'x^, = 0. (3.21) 

(2) (2) 

The gauge invariant variables $ and \1> are the second order scalar perturbations and 

(2) (2) 

z/j and Xij are the second order vector and tensor modes of the metric perturbations, 
respectively. Furthermore, using X"' and the gauge invariant variables for an 
arbitrary tensor field Q can then be defined: 

Wq := (i)g - £xQ,, (3.22) 

:= (2)Q _ 2£x^''>Q - {£y - £1} Qo. (3.23) 

Thus, from Eqs. fl3.22p and fl3.23p . we can see that any variable can be decomposed 
into a gauge invariant and gauge variant part as 

(1) g = (i)Q + i:^go, (3.24) 

(2) Q = (2)Q + 2£x^^^Q + {£y-£]c] Qo. (3.25) 

These decomposition formulae are valid for the perturbations of an arbitrary tensor field 
without knowing detailed information of the background metric gfj^. As a coUorary, any 
equation for the perturbations (e.g. the Einstein equations or the equations of motion 
for the matter fields) is automatically given in a gauge invariant form [231, [2l| due to 
the lower order equations. 

Since we know that any of these equations must be gauge invariant, we only need 
study the gauge invariant parts of Eq. fl3.18l) and not treat the components of lab directly: 
that is, we must only study Cab- However, the choice of the gauge invariant part at both 
first order, Tiab, and second order. Cab, is not unique. This can be seen as follows: 
through the components of Tiab, we can construct a gauge invariant vector field, for 

example, Za ■= —a $ {dr])a + a ^1)1 {dx^)a- Using this we may write 

hab = 'Hab + £zgab — £zgab + £xgab ='■ ^ab + £ xQab, (3.26) 

where we have defined /Cafe := Hab + £zgab and X"' := X"' — Z°' . Thus there are, in 
principle, infinitely many choices of 1-iab, which correspond to the infinitely many choices 
of the gauge fixing. As mentioned in a previous paper [5], the situation at second order 
is more complicated; we simply note here that there exist infinitely many choices of Cab- 
Because of this, it is evident that relating the variables in the MW2009 approach to 
those in the KN2010 approach is more difficult. However, the procedure is the same as 
at first order: compare the components of lab through Eq. (13.181) with Eq. (13. ip . 
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4. Comparison of different formulations 

Having now summarised cosmological perturbations from both the MW2009 approach, 
and that in KN2010, we move to the main point of this paper: linking the two approaches 
and showing the equivalence between them. 



4.1. First order 

In this section we present the comparison between the two approaches at linear order. As 
shown in Eqs. f l2.46p . at linear order, the correspondence between metric perturbation 
variables in the MW2009 and KN2010 approach is quite clear. Here we discuss the 
correspondence of the approaches themselves by choosing two popular gauges: the 
longitudinal (or Poisson) gauge, and the uniform curvature gauge. 



4-1.1. Longitudinal (Poisson) gauge The longitudinal gauge is defined in perturbation 
theory by choosing hypersurfaces of constant shear {au = 0). Thus, in the longitudinal 
gauge, the scalar gauge function au is given by 

au = -cTi = B^-E[, (4.1) 

and this gauge is fully specified (for scalars) by requiring separately that En = (which 
implies that Bu = 0) and is a natural choice. Hence 

Pu = -El . (4.2) 

The remaining two scalar metric perturbations, 0i and ipi, are then given as 

^ = (f)i + H{B, - E[) + (fii - E[y , (4.3) 
^u = ^i-H{B^-E[) . (4.4) 

These are the two gauge invariant Bardeen potentials $ and ^ [T] and, in fact, these 

(1) (1) 

two variables coincide with $ and ^ as defined in KN2010 [5]. 

By including vector perturbations (an extension of the longitudinal gauge generally 
called the Poisson gauge) we must also fix the vector gauge function 71*, which can be 
achieved through the relationship 

lu = j S{dr^ + C{{x=) , (4.5) 

where C\{x^) is an arbitrary constant 3- vector which depends upon the choice of spatial 
coordinates on an initial hypersurface. 

We have thus specified the gauge generating vector, through Eqs. (14. ip . 
(14. 2p . and (14. 5p . In this gauge, the remaining components of the linear order metric 
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perturbation are completely given in the form Eq. (12.441) . Therefore, this gauge-fixing 
corresponds to the choice of the gauge variant part X" in Eq. f l2.22p so that 

X" = pX" = 0, (4.6) 

where V denotes the Poisson (longitudinal) gauge choice. In other words, the gauge 
invariant variables used in the KN2010 approach [5] correspond to the variables 
associated with the longitudinal gauge. 

4- 1.2. Uniform curvature (spatially flat) gauge An alternative gauge choice is the 
uniform curvature, or spatially flat gauge. This amounts to choosing a spatial 
hypersurface on which the metric is jmperturbed by scalar or vector perturbations, 
which requires V^iflat = -^ifiat = and i^ifl^t ~ 0. This gives the gauge transformation 

«lflat = |- , /3lflat = -E, , 7ifiat = • (4.7) 

The remaining scalar metric perturbations are then 

^t=<Pi + i^i+ (jf^ , (4.8) 

B^.t = B^-E[-^, (4.9) 

which are gauge invariant. Thus, in the uniform curvature gauge, the metric 
perturbations as expressed in Eqs. (12. 7p . (12. 8p and (12. 9 p are 

/iflat,7r? = - 2a^</)ifiat , (4-10) 
hfla.ti'q = 2a^ ( -Blflat.i — 'S'lflati ) , (4-11) 



/iflatii = a'^hij . (4.12) 

Now, we compare this with the KN2010 formalism, in which the gauge choice is regarded 
as a choice of the gauge variant part Xa = j^Xa, where J-" denotes the flat gauge choice. 
From Eqs. (^^, (I231D, and Eqs. (^JUi)-(^J^ we obtain 

- 2a' $ +29^(^X^) - 2HijrX^) = -2a^^t , (4.13) 



(1) 

i 

(1) . o (1) 



%\ +D,{j.Xr,) + d^ijrX,) - 2H{jrXi) = 2a'(lhZ,^, - S^,,^, (4.14) 



- 2a" ^ 6ij + a" Xij +A(jXj) + Dj{^Xi) - 2H6ij{^X^) = a^huj (4.15) 

Since the transverse traceless part of both sides of Eq. (14.151) should coincide with one 
another, we may identify 

(1) 

Xij= hiij. (4.16) 
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Further, the trace and the longitudinal part of Eq. f l4.15p is given by 

- +^D'iTX,) - H{j.X,) = 0, (4.17) 

D,{j.X,) + D,{j.Xi)=0, (4.18) 

from which we can choose 

^Xi = 0. (4.19) 

Further, Eq. fHTTD yields 

(1) . . 

jrX, = * • (4.20) 

Then, identifying the divergenceless part of the both sides of Eq. fl4.14p with one another 
gives 

= - ■ (4.21) 



filflat = -i^TXr, = -77 * , (4.22) 



Finally, the scalar parts of Eq. f l4.14p and Eq. fl4.13p yield 

1 1 (1) 

^■^-^'^ - ~H 

' ' (1) / dr,H\ (1) 1 ^ (1) , , 

0ifiat=<f +f ^ +-^5, * . (4.23) 

As expected, this is simply the relationship between the scalar metric perturbations in 
the uniform curvature gauge and the Bardeen potentials, the scalar metric perturbations 
in the longitudinal gauge. 

4-2. Second order 

4-2.1. Longitudinal (Poisson) gauge In order to extend the Poisson gauge to second 
order, we continue as at first order. First, we require that E21 = 0, which fixes the scalar 
part of the spatial gauge as 

(32i = -E2 - ^V-'V-^T' + 7V-'X\ . (4.24) 
4 ' -J 4. 

Then, requiring that i?2£ = sets through Eq. (13.101) . and requiring that the vector 
= can be used to fix the vector part of the spatial gauge transformation, up to 
a constant of integration as at linear order. The second order analogues of the gauge 
invariant Bardeen potentials $ and \& are then given as 

^2 + Ha2i + a2i + au [au + 5Hau + {H' + 2H^) au + 4/701 + 2(f)[] 



+ 2au' {au' + 20i) + ^uk ("i/ + Hau + 20i) + ^uk "i", " 2^1^ - ^ 



k o D ^k 



, I 



-.U 



i^2i = ^2 - Ha2i - + 7 V-'X/ , (4.25) 
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where X^jj is denotes the quadratic first order terms in Eq. f l3.12p using the longitudinal 
gauge transformation components an and ^j^. 

To summarise, the second order metric perturbation has the components 

- 2a^~ 



(Ig)ir] 



a^{'q)i - 2tl'2Sij + h2ij 



(4.26) 
(4.27) 

(4.28) 



As at first order, the choice of the Poisson gauge in the MW2009 approach corresponds 
to the choice 



V 



= 0. 



(4.29) 



in the KN2010 approach. Thus, we again confirmed that the gauge invariant variables 
used in the KN2010 approach [3] correspond to the variables associated with the 
longitudinal gauge, with the correspondence between the variables of the MW2009 
formulation and KN2010 formulation as 



(2) ^ 



(2) ~ 



(2) ^ 
^f-^ ^2, 



(2) ^ 



(4.30) 



4-2.2. Uniform curvature (spatially flat) gauge At second order, the gauge condition 
^|:2 = gives, using Eq. fl3.8p . 

«2fiat = f + ^ [V-^Xt,, - XL.] , (4.31) 

where we get Xflatij from Eq. f l3.12p using the first order gauge generators given above, 
as 



X 



flatij 



+ 4CHj,fc^ifla^t + (4Ciifc + ^lflatj,fc) ^fflatj- + (4Cijfc + ^lfiatj,fc) ^iflat,; 



1 

2 



lfiat(i,j) 



4i/^lflat(jj-)) + 2^lfiat^lflat(i,j)fe + 2^iflatA:,i^lfiat, 



where we define 

^Iflati 



(El,, + Fu) . 



(4.32) 



(4.33) 



Finally, the gauge conditions £'2flat = and -F2fiat ~ ^ enable us to specify the gauge 
functions /32flat and 72fiat completely. Thus, in the uniform curvature gauge, the metric 
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perturbation has the components 



(flat)7?r; 
(flat)-;yi 



2a' 



■'2flat 5 



(4.34) 

(4.35) 

^(flat)ii = a^(??)/i2flatii , (4.36) 

On the other hand, according to the decomposition (13.181) . the second-order metric 
perturbation lab in the flat gauge is given by 

jiab ='■ ^ab + 2£^xjhah + {£ jY — -^^x) 9ah- (4.37) 

Here, the components of jXa are given by Eqs. (I4.19P and (I4.20p . and the components of 
jhab are given by Eqs. fHTTUD - fimi) with the relations fliTTBD . flOTj) . ^TM . and 
Tedious calculations show that the components of the terms 

jE^ah '■= 2£^xj^ab — £'],x9ah (4.38) 



in Eq. fOTj) are 



2a2 



9^ W(l) „ (1),, ^ (1) ^ (1) 

+4H^d^H $^ +4/7^9^/7(^)2 + GHd^H ^ 9^ ^ 

9 ^9 



(4.39) 



17) 



a 



(1) (1) (1) ^ (1) ^ (1) (1) 
8a^i7 ^ A * -3i7 * a^jA * -5ifa^ * /^i * 

(1) (1) (1) (1) 

o (1)(1) n (1) ^ (1) „ (1)(1) 

+AH''dn '^^^Vi , (4.40) 

(1) (1) „ /(i)\2 (1) (1) 

-3D, * L)j ^ -2H\j ( * ) - ^7ii ^ -9,, * 

(i)(i) (i)(i) 
+2HDi ^ z/j +2Mj ^ Ui 

(1) ^ (1) , (i)(i) \ 

+H ^ dr, X^j +2H^ * X.H • (4.41) 



(1) 



(1) (1) 



2a^ 



Comparing different formulations of non-linear cosmological perturbation theory 



20 



Together with Eqs. fl4.34p - p:.36p . we obtain the components of Eq. fl4.37p as follows: 

(2) 



- 2a 02flat 
,2 



2a' $ +2drjTYr, - 2H^Yr, + jE 



mi 



(4.42) 



« ('?) (^-B2flat,i — 'S'2flati 

a^{vi)h 



2 (2) 

a V 



, + D.jY^ - 2HrY, + (4.43) 



(2) (2) 



The trace part of Eq. fl4.44p is given by 







(4.44) 



(4.45) 



and the traceless part is 



(2) 



' — \-/ 2 

o^(^)/^2fla%- = Xij +DirY, + Dj^Yi - -lijO^rYk 



^ii - ^lijT'j^ki, (4.46) 
where jrYj is decomposed as 

^Y, =: D.^Y^L) + tY^v)j, D^tY^v)j = 0. (4.47) 
Eq. K4Q^ is then 

a^(?7)/i2flatji = '^Xij +2 (^DiDj - ^lij^ tY^l) + 2-D(ij-y(y)j) 

+ - ^7ijyjHM, (4.48) 

where we have used the fact that Di is the covariant derivative associated with the fiat 
metric. Taking the divergence of Eq. (14. 48 p . we obtain 



and further taking the divergence of Eq. fl4.49p . gives 



0, 



3 



Substituting Eq. f l4.50p into Eq. fl4.49p . we obtain 

and hence we have 
tY, = 



4: ■' 



(4.49) 



(4.50) 



(4.51) 



(4.52) 
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Eq. (K^ yields 



4 



J^kl 2 ikl I P—'mn 



(4.53) 



T~ _ (2) 

il'2&a.tij — ^ij 

1 

2a 



A-^lkn"''''"jS^mn ■ (4.54) 



On the other hand, the trace part of Eq. fl4.45p gives 



(2) 1 

— qj 

H AH 



A^^D^D^ 



J-^kl e^'^kl'y J-^mn 



Taking the divergence of Eq. f l4.43p and substituting Eqs. fl4.52p . f l4.53p . and fl4.55p . we 
obtain 



2flat 



1 (2) 1 
\[/ 

H 4a2 



3{dr,-2H)A-' + - 



A^^D^D^ 



Qa^H 



(4.56) 



and substituting Eq. f l4.56p into Eq. fl4.43p gives 



(2) 
I/,- 



— {dr, - 2H) [D.A'^D^D^ - 7, ^D^] 



J^kl 2 '~1kl'~i J^mn 



(4.57) 



Finally, through Eq. (14. 42 p . we obtain 

\idr,-H) 



(2) 

^'2flat = $ 



a2 (2) 
— 

H 



(Ji 

1 _ 



:A-^D^D^ 



AH 



J-^kl r^'lkl^ J-^mn 



(4.5^ 



As at linear order, Eqs. f l4.54p . f l4.56p . fl4.57p . and fl4.58p give the relationship between 
between the variables in uniform curvature gauge and the Poisson gauge. 

Before concluding this section, we note the non-uniqueness of the gauge invariant 
variables defined in Eqs. f l2.22p and fl3.18p as discussed at the end of §3.21 In this 
section, we regard the choice of gauge as the specification of the vector fields and 
y", which are the gauge- variant part of the first and second order metric perturbations. 
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respectively. On the other hand, we may also regard gauge-fixing as the specification of 
the gauge invariant parts Hat and Cab of the metric perturbations in the sense discussed 
at the end of §3.21 Actually, we may use jXa, whose components are given by Eqs (14.191) 
and fl4.20l) . as the vector field Za in Eq. f l3.26p . since the components of jrXa are specified 
by the gauge invariant variables. In this case, the first order metric perturbation hat 
is decomposed into gauge invariant part Kab and gauge- variant part X"- as Eq. f l3.26p . 
where ICab and X" are given by 

]Cab = V-ab + £^X9ab. X" = X'^ - j:X\ (4.59) 

Further components of Kab = j^ab are given by Eqs. (I4.10I) - (I4.12I) . JCab is regarded as 
the realisation of the gauge invariant variables for the first-order metric perturbation 
associated with the fiat gauge. 

Furthermore, in the case of the second-order metric perturbation Lb, we define 

J^ab '■= Cab + '^£jrXjhab + {£ jY — £']rx) 9ab (4.60) 

and Eq. (I4.37P is given by 

Lb = Jab + '^£xhab + (y£ y - jY +[X ,X] ~ "^jf j 9ab- (4.61) 

Then, choosing 

the second-order metric perturbation Lb is given by 

Lb = Jab + '2£j^hab + {£y ~ £%) Sab- (4.63) 

The tensor Jab is clearly gauge invariant and Y"" satisfy the gauge transformation rule 
yY^ - xY'' = ^^2) + i^Hh (4.64) 

under the transformation $ = o 3^, i.e., Y"" satisfy the property (13.191) of the 
gauge-variant part of the second order metric perturbation. The components of the 
gauge invariant part Jab are given by Eqs. (I4.34p - (l4.36p and so Jab is regarded as the 
realisation of the gauge invariant variables for the second order metric perturbation 
associated with the fiat gauge. 

At the end of §3.11 we noted that the second order tensor perturbation is gauge 
dependent. Therefore, the expression for the gauge invariant second order tensor 
perturbation will differ depending on the choice of gauge. This is due to the fact 
that in perturbation theory beyond linear order, in which mode-couplings occur due to 
the non-linearity of the system, the notion of the transverse-traceless part of the metric 
perturbation cannot be identified uniquely with the gravitational waves. As shown 
for the second order case in Section 13. Ij the gauge invariant higher order transverse- 
traceless perturbation has contributions from the second order tensor perturbation and 
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other first order metric potentials. In other words, we might say that the second order 
gravitational waves are also generated from the the first order gravitational potentials 
as discussed in literature [2Tj . 

5. Summary and discussions 

In this paper, we have compared and contrasted two different approaches to metric 
based cosmological perturbation theory, and have derived the relationship between 
the standard approach a la Bardeen [1], which was discussed in MW2009 and 
the formalism studied in KN2010 [S]. We started by introducing the basics of 
relativistic perturbation theory, and then presented first- and second-order cosmological 
perturbation theory in the MW2009 and KN2010 approach, respectively. Finally, in §U 
we compared the two formalisms directly for both the longitudinal or Poisson gauge and 
the uniform curvature gauge. 

Both approaches recognise the spurious gauge artefacts present in perturbation 
theory, and adopt different techniques in order to remove them. In the MW2009 
approach, one perturbs the metric tensor and then inspects the gauge transformations 
of the metric perturbations, using these to eliminate the gauge freedom and hence 
construct gauge invariant variables. The KN2010 approach splits the perturbation to 
the metric into a gauge invariant and gauge variant part, and then writes all equations 
in terms of these gauge invariant variables. In the MW2009 approach, the gauge choice 
is made by specifying the gauge generating vector, On the other hand, in the 
KN2010 approach the gauge is defined through the gauge variant vector X" (or, at 
second order, Y""). We showed that the gauge invariant variables that are used in 
KN2010 are equivalent to the usual gauge invariant variables of the Poisson gauge (so 
that the Poisson gauge is specified through X" = = Y"'). When relating KN2010 
formalism to the MW2009 approach in the uniform curvature gauge, we simply obtain 
the usual relationship between the gauge invariant variables in the Poisson gauge and 
those in the uniform curvature gauge. Thus, we have shown that the two approaches 
are equivalent. 

While this result is not necessarily surprising, since both the approaches are 
based upon metric cosmological perturbation theory, showing this equivalence is a 
good consistency check for the KN2010 approach. Furthermore, there may be certain 
problems for which it could be advantageous to use one approach over the other. Having 
shown this equivalence, and having a working knowledge of both theories, may enable 
one to more easily solve the problem in hand. 
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